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ABSTRACT 


A proposed liquid ground state of metallic hydrogen at zero temperature 
is explored and a variational upper bound to the ground state energy is 
calculated. It is shown that the possibility that the metallic hydrogen 
is :i I 'quid around the metastable point (r a 1.64) cannot be ruled out. 

This conclusion crucially hinges on the contribution to the energy arising 
from the third order in the electron* proton interaction which is shown here 
to be more significant in the liquid phase than in crystals. 
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1 . IHTRODUCTION 


An interesting possibility of a zero temperature liquid ground state of 

metallic hydrogen has been recently explored in a calculation which makes 

use of a Jastrow-Slater many particle variational wavefunctlon * to calculate 

the ground state energies of both solid and liquid phases. The symmetric 

part of the wavefunctlon is treated by the Monte-Carlo technique; exchange 

la neglected in the solid and approximated in the liquid by the Wu-Feenberg 
2 3 

expansion It is found that the differences in the energies of the liquid 
and the solid phases varies from 0.1% at r = 1.6 to about 3% at r = 0.8 , 
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(here 4TT/3(r^a^) = 1/n and n is proton or electron density). The solid 

phase seems to be energetically more favorable throughout the entire range 

of densities considered. However, the calculation is based on a model of 

pair-interactions between protons and therefore contains only terms generated 

to second order in the electron-proton interaction. The contribution coming 

from the third order in the electron-proton interaction is known to be slgnl- 
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f leant in the calculation of the band-structure energy * in the solid. In 
view of the small energy difference between the solid and the liquid phases 
it is therefore necessary to estimate the third order term for the liquid as 
well. Furthermore, since in the liquid certain configurations will permit 
three protons to come closer together than they would in a solid, we might also 
expect that the contribution from the term third order in the electron - 
proton interaction may be relatively more important in the liquid phase. 

In this paper we shall first show that a simple one-porometer variational 

O 

wavefunctlon when combined with the Hypernetted Chain (HNC) Integral equation 
can reproduce the energies calculated in Hef. 1 with a 6-parameter variational 
wavefunctlon and the Monte-Carlo technique to within 0.025 - 4.2% and therefore 
provides a very reasonable upperbound. However, precise agreement is not 
necessary in order to provide variational answers to the following questions 
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(a) How much does the third order terra contribute 
to the ground state energy of the liquid? (b) \Vhat are the corrections in the 
liquid state attributable to long wavelength phonons? (c) Is it possible 
to lower the energy o£ the liquid by permitting partial alignment of the 
spins of the protons? 

The calculation described below Is a Judicious combination of variational 
and perturbative methods and is intended to suggest that for certain densities 
the possibility of a liquid metallic phase of hydrogen at zero temperature 
cannot be ruled out. The conclusion hinges on the fact that the third order 
term is significant and is perhaps more so in the livquid. 


2. FORMULATION 

In a sense hydrogen is the simplest metal; its Hamiltonian is known 
exactly: For N protons, N electrons and volume Q we write 
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Here we have denoted the proton coordinates by and the electron coordinates 

by Cr^} . A major simplification takes place® when we realize that there are 
two widely different time scales involved in the problem, allowing us to remove 
electronic degrees of freedom by assuming that at any instant we can consider 
the electrons to be in the ground state corresponding to the instantaneous 
proton configuration. This Born-Oppenheimer adiabatic approximation reformu- 
lates the problem in terras of an effective Hamiltonian of protons. The price 
we pay is that the indirect interaction between the protons, now mediated by 
the electrons, is no longer a simple Coulombic pair interaction but contains 
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many body forces . With electron coordinates now integrated out the total 

s 

Hamiltonian for the protons becomes 


H 
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where E . which is the exact ground statn energy of the interacti ng electrons 
eg ■ ■ 

in a uniform positive background appears as a constant energy, and simply 
drops out of the calculation. In Eq. (2.2) T and V are the parts of the 
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original Hamiltonian of the protons and which are functions of the 

proton coordinates are the electron mediated interactions between protons 
which are generated by adiabatic perturbation theory. Provided Eq. (2.2) 
converges , the procedure is exact within the adiabatic approximation. Most 
importantly, note that to this point we have not made any assumptions regarding 
the positions of the ions; the discussion holds for liquids and crystals 
whether static or dynamic. The precise form of e|^*^^(CrP) cun easily be 


written down 
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S’V(k^)V(-k^) 
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and similarly for the nth order term. Here, 


v(k) ^ 

i k 


(2.4) 


and 
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(2.5) 



4 


Is the exact first order static response of the Interacting electron gas to, 
an external potential. Similarly x ‘ ’ *\+l) exact nth order 

response. In otherwords if we know the nth order response function of the interacting 
electron gas exactly, we would also know exactly these extra many body interactions 
between protons, and we can proceed to diagonalize the proton Hamiltonian. 

The interesting point to note 1s that the rewriting of the original 
Hamiltonian in the form given in Eq. (2.2) splits off a large volume dependent 
terra (order 1 Ry) which does not depend on whether the protons form a 
liquid or a solid and therefore simply drops out of the difference in energies 
between the liquid and the solid phases which is the interesting quantity 
in examining the phase transitions between the two. The uncertainties in 
the electron gas response functions ih„, . - .it ,) will surely affect 

1 2 A'fj. 
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each of the terms (CRj^I) but, once again, they will not Influence too 
greatly the difference in energies. Thus this particular reformulation, Eq. (2.2), 
should be a reliable starting point to calculate the energy difference between 
liquid and solid phases. 

(1) 9 

For X Ot) we shall choose the Hubbnrd-Geldart-Vosko (HGV) form for the 

dielectric function €(k) which is Icno^vn to be of reasonable accuracy at least for 

^3 ^ 2. For X^|(^T we shall make use of the form used by Brovraan, Kagan 

5 

and Holas in which the one body interactions are screened by the HGV dielec- 
tric function. This approximation for x^^^k ,kL,k ) has 

12 3 

teen used extensively and is believed to be reasonably accurate. The 
Hamiltonian can now explicitly be written down^^f we neglect ^or 

n > 4; 
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where, 
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is a large volume dependent term, which is convenient to separate out. In 

(2.7) n is the nunth^r density (NAl) and ^ is the compressibility of the 

uniform interacting electron gas neutralized by a uniform positive background at 

the some density. Note that the terms E^^\[h.]) and V have been combined 

a Sc, pp 

to give 
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an effective linear-response pair potential. Finally the third order term 
is given by, 


( 2 . 2 ) 

K(5i,kj,-Ej-ir3) 

Here X is ; 
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where 0(x) = 1 for x > 0 and zero for x < 0. The remaining parameters are 
given below, 
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(2.13) 


(2.14) 


(2.18) 


(2.16) 


If we take €(k) to be the RPA dielectric functloa then A would precisely be 
the BPA approximation for the three tailed diagram. 

As mentioned earlier the dielectric function e(k) is taken to be of the HGV 
form and la explicitly given as, 


e(Ti) a 1 + 


where 


aF(i1)/7i‘' 


F(“n) a 1 + 


1 - aF(Tl)/(2Tl +g) 
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and “n = k/2I?p, 

Finally, we obtain 


i<j<k ^ 


( 2 , 21 ) 
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wliore Eq is a constant volume dependent term end ue have split o££ the 0^^^ 
term from II ' given by 
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( 2 . 22 ) 
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In lief. H WQ 0 approximated by Eq + H * Wo proceed from this point and 

(25 

shall first attempt to diogonalisse * as well as possible with a one parameter varia- 
tional function which, as we shall see, will give an error of no more than 4% when 
compared tu the calculation of Eef. 1 employing 6 variational parameters. An 

optimum wavefunctlon obtained In this way will be used to calculate the varia- 

(3) 

tlonal bound for the contribution from 0 , 

3. CALCULATIOMAL TECHNIQUE 

In this section we shall outline the method used In calculating the ground 

state energy of the Fermi liquid corresponding to the Hamiltonian given In 

2 3 

Eq. (2.6). A Jastrow-Slater variational wavefunctlon ' 


I (1,2,...N) = D 
— — o 


(3.1) 


will be used to calculate an upperbound to the ground state energy. In Eq. (3.1) 

B 

D Is a Slater determinant made out of plane waves and i]/ is a symmstrlc correlating 
factor designed to take core of the strong Inter-partlcle Interactions. It Is 

2 3 

responsible for a large part of the energy. A subsequent Wu-Feenberg expansion ’ 

then uses an exact transformation to recast the problem into the calculation of 

two distinct parts; Thus we shall set 

B = E- + (3.2) 

j 3 ex 

where E is the exchange contribution and E_ is the eigenvalue of a symmetric 
ground state corresponding to the Hamiltonian. Then 

hcCr^) 


(3.3) 
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v;lierQ ^ In Eq. (3*1) Is chosen to he the eigenfunction of (3.3). The calcula- 
tion of Eq therefore does not involve the antisymmetric factor and results In 
a conslderahly simplified problem. A knowledge of this ^ is then utilized 
to calculate, 


ex 
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which may be calculated by a statistical cluster expansion of the type 


ex 







(3.3) 


where involves n-particle exchange. These terms are easily calculated 

(at least up to the 3rd order) as we shall see below. The entire procedure 

is meaningful when EL is much greater than E and the series in E converges 

o ex ex 

rapidly. We shall see later that the first condition is very well satisfied, 

EL being several orders of magnitude larger than E . However, the second is 
only moderately well satisfied, each term dropping by a factor of 1/3 to 1/S 
of the previous term. 

So far wo have implicitly assumed a paramagnetic ground state, each level 

being doubly occupied in the Slater determinant. However, it is easy to extend 

2 3 12 

the result to a departure from double occupancy ’ ’ . The resulting form for 

E (x) is then 
ex 


Egjj(x) = E^^(x) + E*^^(x) + E°^(x) + ... 


where x is the spin imbalance order parameter defined by. 


N - N 
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Here N (N ) are the numbers of up (down) spins and K Is the total number of 

+ •• 

spins. A non zero value of a will signify a magnetically ordered phase 
Clearly :£ » 1 will represent a ferromagnetlcally ordered phase. Notice that 
does not depend on x. We shall try to determine whether E (x) possesses a 
minimum sT (xm) at a non-zero value of x. It will turn out that the energy 

02S 
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difference AE(x) a E (x=0) - ^ (xm) per particle is small, only ~ 2 x 10 Ry, 
(It is worth noting that this is not small on the scale of a superconducting 
pairing energy.) 

4. VARIATIONAL MimiOD 

From the variational point of view in Eq. (3.2) is conveniently split 
into three parts 




Eg = E^^^ + 4^^ + ^ 

( 2 ) 

The first term, Eg , is calculated by variatlonally optimizing the Hamiltonian 
H^^\[Rjj}) with the many-body Jastrow wavefunction given by, 
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whore , 


u(r) 


A3 


(4.3) 


This wave function is a simplified one-parameter form for that used in Ref. ?. 

The energy functional is minimized with respect to the parameter b at every value 

of r , the resulting wavefunction is then used to calculate the expectation 
a 

value of The obtained in this first order perturbation is also 

a variational bound. The u*(r) expressed in Eq. (4.3) is short ranged and does 
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not Inoludo tho contribution due to tho long wavelength phonons* This Is 

13 

done perturbatlvely with tho help of Choster-Reatto wavofunotlon . The 
relevant formulae are summarized below: 
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where all distances are scaled with respect to the Inverse Fermi wavevector, 

1/kp, including the variational parameter b (b s b^/h^) . In Bq. (4,4), r^ 

denotes the average interparticle distance scaled by the Bohr radius and g^(K) s 

2 3 

g°(r), (r a x/kp) is the pair correlation function defined as; ' 
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Note that ^ is defined in EqB.(4,2)and(4,3). Tho corresponding static structure 

B 2 3 

footor s“(k) is defined by the Fourier transform: ' 


S°(k) a 1 + njd? [g°(r) - 13 
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Finally with tho distance and the wavevector 3> lied, 
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is the screened interaction and e(^'; is the HGV dielectric function. Once 
again all wavevectors are scaled by kp(|kj = shall use the 


u 


2 3 

llypornottQd Chain ApproxiQatlon ' which Is known to be satisfactory for Boso 

2 3 

fluids and has been tested for a variety of Interaction potentials. 

In this approxioatlong^(r) is the solution of the non-linear integral equation 

relating the direct correlation function e<r) to g®(r); 

B 


g®(r) - 1 a c(r) + njdr'cclr-r’ j) [gg(r')-ll, (4.8) 

c(r) a gg(r) - 1 - log ffg(t) + u(r) (4,9) 


The procedure is to solve Bqa.(4.8)and(4.9)for a given value of the variational 

parameter b by a standard numerical procedure and to use the resulting gg(r) 

in Bq. (4.4) to calculate the energy. This process is repeated for a number 

of different values of b to find the optimum go(^)» minimum in 

B 

energy at a given density or r . We then proceed to calculate the contribu- 

s 

tion due to Thus 



-.— 1 —— S^(k.?,-k-q)A(k,T,-k-q) (4.10) 
(q+k) e(q+k) 


where , 


Sg(k,q,-k-q) => 



(4.11) 


and 
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A dlstinot foaturo^^ of tho rosponao function of A(k,q,-k-q) is its singular 
bobavlor when Is -f q a 0: l.o., 

A(k,-k,0) |l^ » kVdl . (4.13) 

This singularity is stronger here than In tho second order response where only 
tho derivative has a logarithmic singularity. This amplification is due to 
the confluence of tho usual second order Kobn anomaly which is always present 
in the third order response and the intrinsic singularity of the third order 
response. It is clear that the integral in Eq. (4.10) can only be defined if 
this fl* '■•ularity la cancelled by other terms present in the integrand. To 
thl . efxect wo prove rigorously in the Appendix the following results 
limjj_^Sg(k, i,“k-£)- ak if limj^_^Sg(k) - ak. Similar results hold when £ 0 and 

Ik+ll - 0. 

Thus It is necessary that S(U) vanish at least linearly with k in the limit 

of small k. Furthermore, any approximation for the three particle structure 

factor must bo such as to preserve this property. One such approximation is 

2 3 

the convolution approximation ’ for the three particle structure factor, an 

14 

approximation that has been extensively tested for soft core potentials and 

14 

in many other situations. Thus we set 

Sjj(k,q,-k-q) - Sg(k)Sg(q) Sg(k+^) ^4.14) 

which clearly has the required property that it vanishes when any of the three 
arguments vanalshes. As is made clear in the appendix this is simply because 
of tho fact that the convolution approximation satisfies all the normalization 
conditions to be required of the probability distribution functions. However, 
as is well known“'^, tho short range wavef unction written down in Bq. (4.3) 
dops not lead to a S_(k) which vanishes as k - 0. This needs to be corrected 
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for the prosenoe, expected physically, cf long range phonons before we can 

evaluote the third order energy given by Eg. (4.10) and (4.14). The procedure 
is 

is almost standard . The Chester and Reatto wavefunctlon is long ranged and 
has the form 
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where we have sealed the distance by kj., l.e. r - x/k^ and x^ is a variational eutoff 

parameter. Here c Is the velocity of sound In this hypothetical Boson system and 
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can be obtained from tho energy, El /N: 
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whore, =j ^ j^nd v^ = (hls^m^). The choice of such a long range wave- 
function leads to a sequence of changes given next. The structure factor 
Sg(k) oal-rulatod with tho short ranged wavefunctlon gets modified to Sj^(k) 
given by 


SgCW 


S^Ck) 


1 * n-s^tk) Oj^jCk) 


(4.17) 


and the corresponding correction in the pair correlation function is 


6g(r) = g°(r) (e"^^^’^- 1) , (4.18) 

where 

gg(r) S5 g°(r) + 6g(r), (4.19) 


and U__(k) is the Fourier transform of U (r) . Finally, 
Lol LR 
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1 + pUj^CU) Sg(k) 


The oorreotlon to the energy is then 
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(4. 21) 


Finally, Eq. (4.10) can be rewritten to obtain the third order contribution 
to the energy. 
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where 9 is the angle between the vectors k and q, 
calculated numerically if S„(q) is known. 
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S (k+q)A(k,q,-k-q) 

e(k+q) 
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Thus Eg ^/N can now be 


5. EXaiMGE CONTRIBUTIONS 


As mentioned earlier the Wu-Feenberg expansion is used to obtain the exchange 
contributions to the energy. total energy per particle is 


E(x)/N = Eg/N + Egjj/N 

= (E^^^ + E^^^ + + E^^ (x)/N (5.1) 

where, E (x)/N is the exchange energy of the Fermions (protons in this case). 
In Eq. (S.l) the energy up to third order in exchange is given by: 


Egjj/N = E^^Ca-,x)/N + Ejg(n,x)/N + Ejg(n ,x)/N + ... 


(5.2) 


where 
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B ^ oJ(l+x)®^®+(l-x)®''®] (5.3) 

1 

^ 2 <n,x)/N a 12Gp{(l+x)®''®J*(y^- | y®+ V) CS(2Is*y)-l]dy 

+ (1-x)®''^ J (y^- i y®+ 4y‘^)CSC2U:y)-l]dy (S.4) 

o 2 F 

and 

* y®23(2;yi2)[S(k;y.3)-lU3(a;yj3)-l)<i7idy2dyJ 

7i<l 

(5.5) 

1/3 

x) and X o (N -N )/N. As mentioned 

+ “ 

earlier our Intention is to compute the ground state energy as a function of 

X. The term E^g is calculated by making the quadratic approximation described in Refs. 

2 and 12 . 

6. RESULTS 

In Fig. 1 we show the dimensionless potential function v^(x) , Eq. (4.7), 

for some typical values of r . In Fig. 2 we show the corresponding pair corre- 

s 

lation functions g„(r). The actual Fermion pair correlation function can be 
6 

2 3 

obtained from these by the Wu-Feenberg expansion ’ , Fermion corrections being 

small in this case. The reason why we have not displayed them is because they 

are not explicitly required in the method of calculating the Wu-Feenberg series 

used here. The structure factor S_<k) corresponding to g„(r) is shown in Fig. 3 

B i3 

for few typical values of r^. It is clear from these plots that there is a 
considerable amount of short range order in liquid metallic hydrogen as compared 
to say liquid helium. One should also note that the interaction potential 
exhibits a strong density dependence. 

( 2 ) 

Table 1 compares our results for , Bq. (4.4), with the calculation 
in Ref. 1. It is clear that our one parameter vax’iational wavefunction gives 


Note that e_ = 
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a reasonably good upperbound. Also shown in the table is tlie 
detailed decomposition of into hlnetlc and potential 

onorglost Wo should oraphaalso that prociso agroemotit between our l->paramotor 

variationol results witii the G-porametor Monte Carlo results, Hof. 1 ,is not 

neoessary since wo are simply interested in an upperbound for the contribution 

arising from tho three body foroos. Those are given in Table 1 along with the 

volume dependent terms . In calculating and we have made use of 

the Nozleres and Pinos interpolation^® formula for the correlation energy of 

electron gas which is consistent with our choice of !IGV dielectric function. 

From Table 1 one can also see that Bq. (*1.21), makes a negligible 

contribution to the total energy. Tho main effect of tho long range phonons 

is to produce on S,.,(k) which vanishes in the limit of small k vvhich, in turn, 

B 

( 3 ) 

allows us to calculate E ^ /N, Bq. 4.22. As noted above the Integral is ill 

conditioned if S,,<k) approaches a non zero value as k goes to zero. 

In Table 2 we have shown, the exchange corrections. It Is seen that a 

partially spin aligned state of pi'otons is in fact favored throughout the 

entire range of densities eonaidored. As mentioned earlier we should be 

P 

cautious about this conclusion since E has been calculated with the help 

(j«3 

2 10 

of the conventional ' quadratic appro.ximatlon, and thus may bo quite 

inaccurate especially for larger values of tho order parameter x. In view 

of the fact that this term is considerably smaller than the rest and that one 

needs a complicated numerical procedure to calculate accurately we have not 

examined it using a more elaborate eomputatioiaal method. We do not believe 
that the x'esults will • change qualitatively. Since the quadratic approximation 

la good in the neighborhood of x » o, the fact that the energy Is lowered for 

non zero values of x can be established although the exact value of x may be 

inaccurate. It is also worth i’emombering that the convergence of Wu-Feenberg 

series is not rigoi'ously established. 

Tlxe total energy for the liquid is compared, Table 3, with the static energies 
the solid phase obtained by Haiamerberg and Ashcroft , .Vote that the static 
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16 

hydrogen could easily be of the order of O.OlRy. Tlie contribution of the 

third order term in the liquid Is more significant than In the solid. For 

example at r » 1.6, the third order energy in the liquid Is -0.0372Ry as 
s 

opposed to -0.0322 calculated by Hommerberg and Ashoroft. The corresponding 

comparison at r » 1.36, yields -0.0326Ry for liquid as opposed to -0.0281 
s 

for the solid Finally, the liquid state energies calculated in this paper 
are a variational upperbound and the exact energy Is expected to be lower. 

Thus one cannot in principle exclude the existence of a liquid ground state 
of metallic hydrogen though it is certainly not established as a preferred 
ground state. 

7. CONCLUSION 

We have investigated the possibility for a liquid ground state of metallic 
hydrogen at zero temperature. We conclude that the possibility of a liquid 
phase near the metastable zero pressure point cannot be ruled out. We have 
found out that the third order terms in the liquid are significantly lower 
than the corresponding ones in the solid and a careful estimate of these terms 
in the solid phase which also incorporates the dynamics of the protons is 
essential to determine the liquid-solid transition (if any). We have also found 
that the contribution to the ground state energy due to the long range phonons 
is negligible though their presence is necessary. An interesting part of our 
calculation is the fact that the energy of this px’oton-electron liquid can bo 
lowered by a partial spin alignment of the protons. 

We would llUe to thank Dr. P. Bhattaoharya and Professor G.V. Chester for 
interesting discussions. This work was supported by NASA, NGR 33-010-188. 



Appendix 


We shall prove that the limiting value of Sg(k,q,-k-q) as any one of the 

wave vector approaches zero from above vanishes provided the static structure 

factor S^Ck) vanishes in the same limit. Strictly speaking this result should 

be considered as a limiting value, defining the function by continuity at 

the origin and true in the thermodsmamic limit. 

2 

Plrat note that , 



S3 


-2+S(k)+S(q)+S(|k+ql)+ 


i foik*? +iq»? -Kk+q)*? ^ ^ ^ 

nJ •*P(rj^,rg,rg) 

drj^drgdrg <A1) 


whore the three particle distribution function P(r^ ,v^ »^3 ^ 


,b2 

N(N-l)(N-2) f»o 

PCri^.rg.rg) g 

n 




(A2) 


Since S_(k,q,-k-q) is invariant with respect to the interchange of its argu- 
B 

ments it is sufficient to prove the result when any one of the wavevectors tend 

» 2 ^ 
to zero, say k - 0 . The following cluster decomposition of P{r^,rg,rg) is 

*4 

exact as long as one does not specify 6P(rj^,rg,rg) j 


P(r^,r2,rg) = n^ Cl+h(r^2)+h(r^3)+h(r2g)+h(r^25^'^^23^-**^^^2^^<^31^ 

+h(rg3^)h(rg2)3 + 8P(r^,?2,?g) (A3) 

where, h(r) “ 62(^5 " 

Then one can easily prove from the normalization of the probability distribu- 

2 

tion functions that 



J‘6P(r^»?3,?g)<i?g « -n®J’h<r^g)hCr23) d?g 


(A4) 


H* 

Now ono oon easily evaluate the right hand side o£ Eq. (Al) for k 0 and 
obtain the stated result. 
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FIGURE CAPTIONS 


Flguro 1 
Flgui'e 2 
Figure 3 


Table 1 

Table 2 
Table 3 


v^(r) for some tynleal valutas of r 
o 8 

g-(r) for some typical values of r 

o 8 

SqCR) for aome typical values of r^ 

TABLE CAPTIONS 

Boson part, E^, of the ground state energy. Eg (MC) is 
the Monte-Carlo results of Ref. 1. All energies are e.spres^ed 
in units of Rydbergs. 

Exchange contribution to the ground state energy. All energies 
are expressed in units of Rydbergs. 

Comparisons of the ground state energies of the liquid (E(x)/N) 
and the solid phases (E (RA)/N: Hammerberg and Ashcroft, Ref. 4). 

All energies are expressed in units of Rydbergs. SC: Simple 

cubic*, BCC: Body centered cubic; FCC: Face centered cubic. 
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TABLE 2 


V 


8 


K 


E^^(S)/N 


0.50 

0.580 

0.00263 

o.eo 

0.579 

0.00102 

1.20 

0.582 

0.00045 

1.30 

0.SB5 

0.00039 

1.36 

0.587 

0.00035 

1.40 

0.588 

0.00033 

1.45 

0.561 

0.00031 

1.50 

0.903 

0.00029 

1.55 

0.S6S 

0.00027 

i.eo 

0.508 

0.00026 

1.70 

0.603 

0.00023 

1.80 

0.607 

0.00021 


TABLE 2 




E®(HA)/N 



SC 

Fce 

BGC 

0.50 




0.30 




1.00 

-0.71188 

-0.71929 

-0.71819 

1.20 

-0.93796 

-0.94019 

-0.93902 

1.25 

-0.96842 

-0.98961 

-0.96843 

1.30 

-0.90217 

-0.99242 

-0.99122 

1.36 




1.50 

-1.04104 

-1.03818 

-1.03693 

1.60 

-1.04739 

-1.04345 

-1.04222 

1.65 

-1.04803 

-1.04338 

-1.04209 


1.70 

1.80 


E<s)/N 

3.36309 

• 0.08811 

- 0.91901 

- 1.00150 

-1.03385 

- 1.04322 

- 1.04509 

- 1.04178 











